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The black hole physics techniques and results are applied to find the set of the ex- 
act solutions of the SU(3)-Yang-Mills equations in Minkowski spacetime in the Lorentz 
gauge. All the solutions contain only the Coulomb-like or linear in r components of 
SU(3)-connection. This allows one to obtain some possible exact and approximate solu- 
tions of the corresponding Dirac equation that can describe the relativistic bound states. 
Possible application to the relativistic models of mesons is also outlined. 

1. Introductory Remarks 

The present paper is motivated by the part of hadronic physics studying mesons. 
As is known, mesons within the hadron physics are the one of central sources of 
information about the quark interaction. Since nearly all known mesons are bound 
states of a quark q and an antiquark q' (where the flavours of q and q' may be 
different) then actually all modern meson spectroscopy is based on one or another 
quark model of mesons. Referring for more details, e. g., to the recent up-to-date 
reviewel and references therein, it should be noted here that at present some gen- 
erally accepted relativistic model of mesons is absent. The description of mesons 
in quark models is actually implemented by nonrelativistic manner (for example, 
on the basis of the Schrodinger equation) and then one tries to include relativistic 
corrections in one or another way. Such an inclusion is not single- valued and varies 
in dependence of the point of view for different authors (see, e. g. RefsB and ref- 
erences therein). It would be more consistent, to our mind, building a primordially 
relativistic model so that one can then pass on to the nonrelativistic one by the 
standard limiting transition and, as a result, to estimate the relativistic effects in 
self-consistent way. 

As follows from the main principles of quantum chromodynamics (QCD), the 
appropriate relativistic models for description of mesons as relativistic bound states 
of quarks should consist in considering the solutions of Dirac equation in a SU(3)- 
Yang-Mills field representing gluonic field. The latter should be the solution of the 
corresponding Yang-Mills equations and should model the quark confinement. Such 
solutions are usually supposed to contain at least the components of the mentioned 
SU(3)-field which are Coulomb-like or linear in r, the distance between quarks, and 
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in what follows we call these solutions the confining ones. No general recipe of 
obtaining such solutions exists at present. In the given paper for this aim we shall 
employ the techniques used in Rcfs.13 for finding the U(N)-monopole solutions in 
black hole physics and, to analyse the corresponding Dirac equation, we shall use 
the results of RefsJil about spinor fields on black holes. 

Further we shall deal with the metric of the flat Minkowski spacetime M that we 
write down (using the ordinary set of local spherical coordinates r, <p for spatial 
part) in the form 

ds 2 = g^dx* <g> dx v = dt 2 - dr 2 - r 2 (d$ 2 + sin 2 tidip 2 ) . (1) 

As is known, such a metric can be obtained from the Schwarzschild black hole metric 
when the black hole mass is equal to 0. Besides we have \8\ = \ det(g^)\ — (r 2 sin d) 2 
and < r < oo, < < 7r, < tp <2tt. 

Throughout the paper we employ the system of units with Ti = c = G = 1. 
Finally, we shall denote L2{F) the set of the modulo square integrable complex 
functions on any manifold F furnished with an integration measure while L^F) 
will be the n-fold direct product of L2(F) endowed with the obvious scalar product. 

2. Dirac equation 

To formulate the Dirac equation needed to us, let us notice that the relativistic 
wave function of meson can be chosen in the form 

with the four-dimensional spinors ipj representing j-th colour component of meson. 
Under this situation, if denoting S(M) and £, respectively, the standard spinor 
bundle and three-dimensional vector one (equipped with a SU(3)-connection) over 
Minkowski spacetime, we can construct tensorial product S(AI)<S>£,- It is clear that 
ip is just a section of the latter bundle. Then the corresponding Dirac equation for 
ijj can be obtained with using the construction of Dirac operator T> with coefficients 
in £ and may look as follows 

Zty = MoV>, (2) 
where /io is the reduced relativistic mass which is equal, e. g., for quarkonia, to 
one half a constituent mass of quarks forming quarkonium, but definition of /io 
is, generally speaking, notsingle-valued and requires the specification within the 
concrete problem (see Ref.Ej). Under this sitiation the coordinate r makes sense of 
the distance between quarks. 

From general considerationsQ'ErQ the explicit form of the operator V in local 
coordinates x 11 on a 2fc-dimensional (pseudo)riemannian manifold can be written as 
follows 

V = i( 7 e ® / 3 )££(^ ® h - \unabl a l b ®h- igA„), a<b, (3) 
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where A = A^dx^, A^ = A^T C is a SU(3)-connection in the bundle £, J3 is the 
unit matrix 3x3, the matrices T c form a basis of the Lie algebra of SU(3) in 
3-dimensional space (we consider T a hermitean which is acceptable in physics), 
c = 1,...,8, <g) here means tensorial product of matrices, g is a gauge coupling 
constant. Further, the forms uj a b = ^>^abdx^ obey the Cartan structure equations 
de a = uj a b A e b with exterior derivative d, while the orthonormal basis e a — e®dx^ 
in cotangent bundle and dual basis E a — E^d^ in tangent bundle are connected by 
the relations e a (Eb) = 6%. At last, matrices 7° represent the Clifford algebra of the 
corresponding quadratic form in C 2 . Below we shall deal only with 4D lorcntzian 
case (quadratic form Q13 ). For this we take the following choice 

for 7° 

where a b denote the ordinary Pauli matrices. It should be noted that, in lorcntzian 
case, Greek indices fi, v, ... are raised and lowered with g^ u of (1) or its inverse g^ v 
and Latin indices a,b, ... are raised and lowered by r\ a h = if b = diag(l, -1,-1,-1), so 
that effigy = if\ E£E»g^ = Vab and so on. 

Using the fact that all bundles over Minkowski spacctime are trivial and, as a 
result, they can be trivialized over the chart of local coordinates (t, r, tp) covering 
almost the whole Minkowski manifold, we can concretize the Dirac equation (2) 
on the given chart for ip in the case of metric (1). Namely, we can put e° = dt, 
e 1 = dr. e 2 = rdfl, e 3 = rsinfldip and, accordingly, E = d t , E\ = d r , E 2 = d$/r, 
E3 = dp/^siwd). This entails 

u>i2 = — di9,o;i3 = — sin $ dip, lu 23 = — cosfldip. (5) 

As for the connection A^ in bundle £ then the suitable one should be the confining 
solution of the Yang-Mills equations 

dF = F A A - A A F , (6) 

d*F = *F/\A — A/\*F+J (7) 

with the exterior differential d = dtdt + d r dr + d^d'd + d v dip in coordinates t, r, ip 
while the curvature matrix (field strentgh) for ^-bundle is F — dA + A A A and * 
means the Hodge star operator conforming to metric (1), J is a source. It is clear 
that (6) is identically satisfied — this is just the Bianchi identity holding true for 
any connection (see, e. g., Refs.i) so that it is necessary to solve only the equations 
(7). Introducing the Hodge star operator * on 2-forms F — F^ u T a dx^ A dx v with 
the values in the Lie algebra of SU(3) by the relation (see, e. g., Ref.i) 

(F^dx* Adx")^F c ^dx a f\dx (3 ) = (g^ a g v0 -g^g va )F^F^y/5dx°A- ■ -Adx 3 (8) 

written in local coordinates x M [there is no summation over a in (8)] , in coordinates 
t,r,d,p we have the relations 

*(dt A dr) = —r 2 sin$d$ A dip , *(dt A dd) = sinddr A dp , 
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*{dt A dtp) = -^-^dr A d-d , *{dr A d'ff) = sin ddt A dtp , 
sin v 

*(dr A dip) = — ^—dt A d-d , A dp) = „ 1 n cfr A dr , (9) 
sin u r z sin W 

so that * 2 = ** = — 1, as should be for the manifolds with lorentzian signature.0 

Besides the sought solutions are usually believed to obey an additional gauge 
condition. In the present paper we take the Lorentz gauge condition that can be 
written in the form div(j4) = 0, where the divergence of theJ^ie algebra valued 
1-form A = A c ^T c dx^ is defined by the relation (see, e. g. Refs.Q) 



div(A) = -=d^\6\g^A„) . (10) 
V\ s \ 



3. Confining solutions 

Now we can use the techniques of RefsJi to find a set of the confining solutions 
of Eq. (7). The essence of those techiniques consists in systematic usage of the 
Hodge star operator. Let firstly J = and let us put T c — A^, where A c are the 
Gell-Mann matrices (whose explicit form can be found in Refs.0), and in detail let 
us write out the addend A^ = A^T C of the operator V of (3) 



A^K 



J_/l8 

A], + iAl 



4 + iA l 



A 1 

-A 3 - 
A 6 

A* 



iAl 
iAl 




(11) 



Then it is naturally to put A c ^ — with c = 1,2, 4, 5, 6, 7 and to obtain from (2) the 
following system of Dirac equations for colour components ipj 



<9 M - -w M afc7 a 7 



if) 



4 + "7=4 



(12) 



2^Vh7 7 



'.9 



4 + 77=4 



02 = MO ^2 , 



(13) 



<9 M - 2 W ^fc7°7 b ~ *5 



v/3 " 



■03 = MoV>3 



(14) 



Further for to avoid unnecessary complications with the Lorentz condition, it is 
simpler of all to put A \ = 0. After this we search for the solution of (7) in the 
form A = A t {r)dt + A v (r)dip with A t , v = A 3 V X 3 + Af^Xs- It is then easy to check 
that A A A = 0, F = dA = —d r A t dt A dr + d r A v dr A dip and, with the help of the 
relations (9), *F — r 2 smdd r A t d-d A dp — -^X^d r A v dt A d-d. From here it follows 
that *F A A — A A *F and Eq. (7) will be equivalent to the equation d * F = 
while the latter yields 



d r (r 2 d r A t ) = 0, d 2 A v = , 



(15) 
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and we write down the solutions of (15) in the combinations that are just needed 
to insert into (12)-(14) 



2 as a 2 

i 

1 



-i§*-T +A *> (16) 



A% + -±=A B V = hr + Bi 
-A% + ±=A% = -(h + b 2 )r - (B, + B 2 ) , 

-l=A% = b 2 r + B 2 (17) 

with some constants aj , Aj , bj , Bj parametrizing solutions. 

Another class of the confining solutions of (7) can be obtained with the aid of 
the ansatz A = A t (r)dt + A v (r, $)dtp, i. e. now we consider the component A v 
depending also on §. Again we shall have A /\ A = 0, F = dA = —d r A t dt A dr + 
d r A v dr A dip + d$A v dd A dip and *F = r 2 sini!)d r Atd-d A dip - ^^d r A v dt A d-d + 
Hiin~g dfl^4yCfc A dr. Then anew *F A A = A A *F and Eq. (7) is converted into 
d * F = that entails 

d r (r 2 d r A t ) = 0, (18) 

r 2 d r A v + sin^ f -^— a d#A v ) = . (19) 
\ sin v J 

It is clear that Eq. (18) gives the same solutions of (16). We shall not here discuss 
the general form of the solution for Eq. (19) and only write out the possible solution 
of it which is useful to us in the present paper in the form slightly modifying (17) 



A% + -^=A% = -K x costf + hr + B x 



-A% + —A 8 = {Kx + K 2 ) cost) - (6i + b 2 )r - {B x + B 2 ) , 
V3 

- A A% = -K 2 cos # + b 2 r + B 2 (20) 

with some real constants Kj, bj, Bj parametrizing solution. 

Finally, we shall adduce one solution of (7) with a source. It is given by relations 
of (16) for At while for A v we have 

A v = I 3 [(br + B)sintf + Kcosi?] (21) 

with vectors b = {h,-h - 62,62), B = {B 1 ,-B 1 - B 2 ,B 2 ), K = {-K^Kx + 
K 2 , —K 2 ). The conforming source is J = — p g I3 (br + TS)dt A dr A d$ and since 
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*J = j= j^dx^ = ) v dcp with j v = ~ ri \ nrd I 3 (hr + B), div(j) = 0, we can interpret 
j as the colour current giving rise to the corresponding linear interaction. 

At last, it is easy to check that all the solutions obtained satisfy the Lorentz 
gauge condition. 



4. Spectrum of bound states in the Coulomb-like case 

Returning to the Eqs. (12)-(14), we should note that when inserting the confin- 
ing solutions gained into them the variables r and i? cannot be, generally speaking, 
separated. Let us take, for example, the solutions of (16) and (21) for the Yang- 
Mills equations (7) with the source J described late in Sec. 3 (where we at first put 
Kj = 0) and for all three equations let us employ the ansatz 



F 32 {r)a^ 3 {d^) 



,j = 1,2,3 



(22) 



with a 2D spinor $j = I ^ j. Then, in the way analogous to that of Refs.Q, we 
can get, e. g., from (12) the system 

1\ 1 



l 

r 

with ci = u>i 



-V - a 2 g h + 



-Fn$i = i(/x - ci)-Fi2$i, 
r r 



1 



-V - a 2 g b l 



Bi 

r \ r 

g{—a\/r + A\) while 

1 



1 



-F 12 eri$i 



-iiVo + cJ-FuCT!®! (23) 
r 



. a 0-i<7 3 d v - -CT 2 (T3COSl9 

smw V 2 



(24) 



is the euclidean Dirac operator on the unit sphere S 2 . It is not complicated to 
check that at b\ ^ 0, B\ ^ the variables r and -d are not separated. Under this 
situation we shall at first restrict ourselves to the case bj = Bj = since under the 
circumstances we can use the results of Rcfs.S for to solve Eqs. (12)-(14) exactly. 
Primarily let us put also Kj — and then we shall modify results to take into 
account the case Kj ^ 0. For all three equations we employ the ansatz (22) and in 
the way analogous to that of Refs.Q, we can obtain from (12)-(14) the systems 



d r + — Fji = iifio - Cj)F. 



:;2- 



r 



F j2 = -i(fi + Cj)Fj! 



(25) 



with an eigenvalue \j for the eigenspinor $j of the above operator 2?o, Xj = ±(lj 
1) G ^\{0} , lj = 0, 1, 2... (for more details see Refs.l). Besides 



C2 = <JJ 2 



ai + a 2 



(A 1 + A 2 ) 
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c 3 = ^3-.g(-^ + ^ 2 ) , (26) 

so that the energy spectrum e of meson is given by the relation e = u>i A- u>2 + ^3- 
The explicit form of $j is not needed here and can be found in Refs.a For the 
purpose of the present paper it is sufficient to know that spinors <£>,, can be subject 
to the normalization condition 

| 2 + \9j2\ 2 )wa.#d0dtp = 1 , (27) 



i. e., they form an orthonormal basis in L|(S 2 ). Further we restrict ourselves to 
the case j = 1 because for j =2,3 the considerations are the same. 

Let us employ the ansatz Fn = \J fio — (wi — gAi)r ai e~ l3ir [fu(ri) + /i2(ri)], 
F u = Vmo + (^i - glQr^e-^ [/ n (n) - / ia (n)] with ax = ^X{-g^a{, fa = 
\J /j,q — (uii — gAi) 2 , r\ — 2/3\r. Then, inserting the ansatz into (25), adding and 
subtracting equations give rise to 

/Wil+*l/ll + -Zl/l2 = 0, 

/W12 - Pmfu + Yifxi + Z 2 f u = , (28) 

where prime signifies the differentiation with respect to r\, Y\2 = <^iPi T 9 a i ( w i — 
gAi), = Ai/?i ± gai[j,Q. From (28) one yields the second order equations in r\ 

r 1 /{' 1 + (l + 2a 1 -r 1 )/{ 1 -^/ u =0, 

Pi 

n/(' 2 + (l + 2«! - n)/[ a - fl + ^) /i2 = , (29) 

that are the Kummer equations (see, e. g. Ref.i) and their only solution finite 
at and at infinity not strongly increasing is the Laguerre polynomial If m (r\) 

(that always may be normalized to be equal to 1 at n = 0) with p = 2a% and 
-nW = Y 1 /(3 1 = 0, -1, -2, ... If putting 

/n = C n L p (1) (r 1 ),/ 12 =C 12 ^ (1) (n) 

then at n = from Eqs. (28) it follows CnYi + C 12 Zi = 0, C n Z 2 + C 12 Y 2 = 
and since Y{Yi — Z\Z 2 =0we obtain C\2 = — Cn^- = —C\\^-. For to describe 
relativistic bound states we should require ip E L 2 2 (M. 3 ) at any ( e 1 and one can 
accept the normalization condition for Fj±, Fj 2 in the form 

f°° 1 
J (\F j i\ 2 + \F j2 \ 2 )dr = - (30) 

with taking into account the condition (27) so that C\\ can be determined from 
the relation (30). The condition Y\j{3\ — —n^ 1 = signifies Y\ = that entails 
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Z\Z 2 = 0. Then at Ai > 0, Z\ = 2gai[io > (we consider a% > 0), Z2 = 0, 
Cxi = and n[ 1} = is admissible, i. e., n { 7 1] = 0, 1,2.... At A x < 0, Yi = Z x = 
but Y 2 = 2501(071 — gAi) ^ 0, Z 2 = —2gai^o 7^ 0, C12 7^ and 717- = is not 
allowed, i. e., n£ = 1)2, .... Under the situation the equality Y\/ [3\ — —rir gives 
the spectrum (provided that w\ > gA\ > 0) 



0J1 = gAi + /i 



(n?> + y/Xf - g 2 aiy 



-1/2 



(31) 



Further, if considering a 2 > 0,w 3 > gA 2 > 0, cj 2 < with \uj 2 \ > g{A\ + A 2 ) then 
acting in line analogous with the above we have 



^3 = .9^2 + Mo 



2 2 

9 ^ 



1 -1/2 



(n (3) + v ^|3^|)2j 



u 2 



Mo 



1 



g 2 {ai +a 2 ) 2 



{n? ] + ^A 2 - 3 2 (ai+a 2 ) 2 ) 2 



-1/2 



and the spectrum of meson is obtained in the form 

'1/2 



£ 

Mo 



1 + 



S 2 «i 



1 



g 2 (ai +a 2 ) 5 



(32) 



l -1/2 



2 2 

5 a 2 



(n?> + ^\ 2 2 -g2(a 1+ a 2 r)2\ 
-1/2 



(33) 



where the numbers 



,(2,3) 



are subject to the same conditions as for n r 



(i) 



5. Influence of the Dirac-like monopole configurations of gluonic field 

If Kj ^ then according to the results of Refsfl we shall get the similar re- 
lations providing that Kj = kj/g with kj £ Z, i. e., kj are integer numbers. 
But now we should consider the spinor of (22) to be the eigenspinor <5>j of 
the twisted euclidean Dirac operator T>k on the unit sphere S 2 , respectively, with 
the Chern numbers k — k±,— (ki + k 2 ),k 2 and the eigenvalues Xj should be, ac- 
cordingly, replace d by Ai = ±y/(h + l) 2 - k'{, h > \ki\, A3 = + I) 2 — |o, 
^3 > \k 2 \, A 2 = ±y/(l 2 + l) 2 — (k~i + k 2 ) 2 , l 2 > |fci + fc 2 | (for more details see Refsia). 
Physically the corresponding configurations of gluonic field describe the Dirac-like 
monopole ones with magnetic charges, conformably, P\ — k\/g, P 2 = —(ki + k 2 )/g, 
P3 = k 2 /g but the total (nonabelian) magnetic charge of the given configurations 
remains equal to Pi + P 2 + P3 = 0. A number of authors (in particular, those 
who develop lattice approach) insist on important role of such configurations in 
a possible mechanism of the quark confinement (see, e. g. Refs.u and references 
therein). 
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6. Spectrum of bound states in the Coulomb-linear case 



We can get the approximate solutions of Eqs. (12)-(14) if considering a 2 $j = $j 
for the eigenspinor &j of (22). As follows from explicit form of those spinors found 
in Refs.u this can be fulfilled only approximately and it is clear that when doing 
so we make an error retaining eigenvalues Xj of the euclidean Dirac operator T>o 
on the unit sphere § 2 instead of the eigenvalues of a less symmetric operator on 
S 2 whose form is unknown explicitly. As will be seen below, however, this may 
seemingly be compensated by the choice of parameters Bj of models so that the 
solutions obtained further may be considered as the almost exact ones. Let us 
employ the solutions of (16) and (21) for the Yang-Mills equations (7) with the 
source J described late in Sec. 3 where at first we put Kj — as well. Using the 
given solutions and the ansatz (22) we shall, in line with the above, get, e. g. from 
(12), the system (23) which (having accepted 02$., « $j) yields 



fl a. Al ft, _l Bl 

Or H g[h-\ 

r \ r 



O r h 9 [bi H 

r \ r 



Fu = «(Mo - ci)Fi2, 



Fn = -i(jM)+ci)Fi 1 . (34) 

Now we employ the ansatz F n = Ar^e-^lfuir^+fuin)], F 12 = iBr ai e-^ r [fn{ri) 
fa(ri)] with ax = ^(X 1 -gB 1 ) 2 -g 2 al, pi = - fa - gAtf + g 2 ti(, A = 
gb\ + B = ft a + uji ~ gAi, r\ — 2(3\r. After this, inserting the ansatz into (34), 
adding and subtracting equations entail 



nABf' n + Yi/n + Zi/ 12 = , 
nABf[ 2 - nABfa + K2/12 + \Z 2 - j^ ABr i ) /a = 



(35') 
(35") 



where prime signifies the differentiation with respect to r±, Y\. 2 — [cxiPi Tffii( w i 
gA 1 )+ga 1 b 1 ]B±g 2 a 1 b 1 A, Z 1<2 = [(Ai - gB ± )A ± ga U i )]B ± g 2 a 1 b 1 A and 
ZiZ 2 = 0. 

From (35'-35") one yields the second order equations in n 



ri /(' 1 + (l+2a 1 -ri)/{ 1 +n( 1 )/ : 







n/r 2 + 



2a t -ri)f{ 



a — br\ 

with a = Z 2 , b = gbiAB/fii, k = AB/Y± and 

ghZi - fcYi 



a — br\ 



l /i 



,(D 



PxAB 



(36) 
(37) 

(38) 



It is clear that according to (36) we should choose / n ~ L p (1) (ri) with the Laguerre 
polynomial V {1) {r\) (p = 2a±) if rir = 0, 1,2... 
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A little bit below we shall show that one can choose ./Wn) ~ P P a)( r i) with 

some polynomial P p (1) (ri), consequently, we shall gain that ipi £ L 2 (M. 3 ) at any 

t £ R and, as a result, we obtain from (38) that uii should be determined from the 
equation 

[g 2 a\ + (nW + ai) 2 ](^i - gAi) 2 + 2(Ai - s-Bi)s 2 ai&i(a;i - gAi)+ 

[(Ai - .gBi) 2 - («W + ai) 2 ]g 2 b\ - »l(nW + ai f = , (39) 

that yields 

A {Xi - gB 1 )g 2 a 1 b 1 - VXT 

w i = 5^1 — —(I) ^ ( 4 °) 

g 2 a 2 + (n| ; + ai) 2 

with Xj = (Ai -gB x ) 2 g^a\b\ - [g 2 a\ + (n, (1) + a- 1 ) 2 ]{[(A 1 -gB x ) 2 - (n^ + ai ) 2 ]g 2 b 2 - 
^{rir + «i) 2 } and the sign of \[X\ is chosen from the fact that the expression 
(40) should pass on to (31) at b\ = B\ = 0. The similar relations will also hold 
true for u>2,3 with replacing ai,Ai,bi,Bi — > a 2 ,A 2 ,b 2 ,B 2 respectively for W3 and 
a±,Ai,bi,Bi — ► — (01 + a 2 ), — (Ai + A 2 ), -(61 + b 2 ), —{B\ + B 2 ) respectively for uj 2 
and the spectrum of meson will be 

2, 



(A x - gB x )<?a x b x - VX 



g 2 a\ + (n^ + «i) 2 
[A 2 + g{B x + B 2 )}g 2 ( ai + a 2 )(h + b 2 ) + y/Xi (A3 - gB 2 )g 2 a 2 b 2 - ^X~ 2 



, (41) 



5 2 (ai + « 2 ) 2 + («? } + « 2 ) 2 g 2 a\ + (n, (3) + a 3 ) 2 

where 012.3 and ^2,3 are obtained from ati, X\ by the same above replacements while 
the signs of \fXj are chosen to correspond to the expression (33) in the Coulomb- 
like case bj = Bj = 0. We can see that all dependence of Xj is expressed through 
the combinations Xj ± gBj and the approximation error mentioned early in this 
section may probably be compensated by choosing Bj in a proper way. 

Now let us describe the above polynomial P p (1) (ri) which should be the solution 
of the equation (37). Instead of solving (37), however, we use the fact that the sought 
polynomial should also obey the first order linear equation (35") with /n(ri) = 
L p (1) (ri) which allows us to write down P p (1) (ri) in the form (putting = n, 
T\ = x for brevity) 

PZM=x-°f£r(%-^)e^t)dt (42) 

with a = 2a>i + n — gb\Zxj{fi\AS). Further we can employ the relationS 
r x x a (l + o) 

/ t^e^L^dt = — , 2*2(0!, 1 + p + n;a + 1,1 + p;-x) ,Rea,x > 



with the generalized hypergeometric function 2*2 of order (2,2) while the Pohgam- 
mer symbol (a)k = a(a + l)...(a + k — 1) = T(a + k)/T(a) with the Euler gamma 
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function T(z), the rclatiorfl 

2 F 2 (a,l + p + n;a + l,l + p;-x) = E ttt^T C n T^TT\ ®(a + k,a + k + l;-x) 



k=0 



(l+ P ) k n (a+l) k 



with the Kummer function $ = 1F1 and the relational <$>(z, z+1; —x) = zx z "f{z, x) 
with the incomplete gamma function 7(2, x) to come to the expression 



P P(x) = il + p)n x ~ a e x 



E 



j^-l( a + k+l,x)- ^7(" + k, x) 



(43) 

Now, using the functional equation 7(2+1) = 27(2)— x 2 e _x , it is not complicated 
to see that 7(0; + k + 1, x) = (a + k)j(a + k)~- x a+k e~ x , j(a + k,x) = (a)kj(a, x) — 
Ei=i( a + i)k~iX l ~ 1 ]x a e~ x and then one can show with the help of (38) and of 
equality Y{¥ 2 - Z X Z 2 = that 



(a + k) ; — 

pV AB 



{a) k = . 



(44) 



At last, employing the relatior£3 

(l+/0n V-(-l)***w. 



E 



we come to the final form of the sought polynomial 

Pi pi 



TO 



(l + p)„^ (-1) 







E^ a+ *) fe - 



^(^^ = 1,2,3. 



(45) 



Accordingly, we can now put fn — CuL^x), / 12 = C\ 2 PP(x) and from the system 
(35'-35") at n = x = it follows the relation between the constants Cn, C12 



C 



11 



ZiP£(0) _ Y 2 PP(0) 

'^12 — ^ To77vT^12 • 



n i£(o) 



Z 2 L£(0) 



(46) 



For completeness let us adduce the normalized radial part of ipi -component for the 
meson wave function of ground state (where n = and then C\\ = C\ 2 = C) 



F 11 =Cr^e-^ r [l-^j (<?&!+ p\), 
F 12 = iCr ai e~^ r [ 1 + 9 -^j ( Mo + wi - ffAi) 



(47) 
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where C is determined from the relation (30) with the help of formulaEI J °° x a 1 e px dx 
T(a)p~ a , Re a,p> 0, which entails 



C 2 



i - ^) 2 (gh + x f + (i + I 1 ) 2 (Mo + - .g^i) 2 



r(2ai+l)(2/3 1 )- (2Ql+1) 



(48) 

It is evident that the similar expressions will hold true for the corresponding ra- 
dial parts of ■(/'2,3-components with taking into account the replacements described 
above. 

Finally, it should be noted that the influence of the Dirac-like monopole configu- 
rations for gluonic field (when Kj ^ 0) can be treated by the same manner as in Sec. 
5 if taking 02 $j ~ &j for the eigenspinor of the twisted euclidean Dirac operator 
T>k on the unit sphere § 2 with the conforming Chern numbers k = kj_, — (fci+fe), ki. 



7. Concluding remarks 

So we have seen that the black hole physics methods and results may really 
help to build relativistic models of mesons that are parametrized by constants 
fj,0i 5S a ji A 3 bj ;, Bj, Kj which should evidently be determined by comparing to ex- 
perimental data. Even in the case of pure Coulomb interaction these models give 
essentially different spectrum of bound states in comparison with the naive picture 
of threefold positronium-like spectrum accepted, e. g., in modern quarkonium spec- 
troscopy. It seems to us that the exact equations derived in the paper may serve as 
a basis to construct miscellaneous approximate approaches in the spirit of that of 
Sec. 6. We hope to continue further study of the questions raised here. 
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